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The problem of large-scale correlations of particles produced in high-energy col- 
lisions is discussed. Among them are, e.g., those correlations which lead to 
ring-like and elliptic flow shapes of individual high-multiplicity events in the po- 
lar+azimuthal angles plane. The wavelet method of analysis which allows to dis- 
entangle various patterns is proposed and applied to central Pb-Pb collisions at 
158AGeV. 

Various correlations of particles produced in high-energy collisions are 
known. Especially well studied are the two-particle correlations due to the 
decay of two-particle resonances and due to Bose-Einstein effect for identi- 
cal particles. The common correlation function technique is well suited for 
this purpose. It is more difficult to apply it if several particles are involved. 
Nevertheless, correlations leading to clusters (or mini-jets, non-reducible to 
resonances) and to jets have been studied as well, however mostly in e + e~- 
collisions where jet and subjet structures are clearly visible. Otherwise one 
has to apply the averaging procedure to get knowledge of dynamical correla- 
tions from studies of the moments of various distributions and their behavior 
in different regions of the phase space (for recent reviews, see Em). It has lead 
to understanding of some global features of QCD Bu. 

However, the genuine multiparticle correlations originating from some 
collective effects may appear in individual events. There exists a method 
of the wavelet analysis which allows to recognize patterns due to correlations 
at different scales and damp down the statistical noise even if it is very large 
in an initial sample. 

Before delving into applications of this method, let us consider examples 
of possible correlations which one would seek for. Everybody knows about 
Cherenkov radiation of photons by a bunch of electrons traversing a medium 
whose refractivity index exceeds 1. These photons form a ring in the plane 
perpendicular to electrons motion, i.e. they are emitted at a definite polar 
angle. As a hadronic analogue, one may treat an impinging nucleus as a bunch 
of confined quarks each of which can emit gluons when traversing a target nu- 
cleus. I speculated □ about possible Cherenkov gluons relying on experimental 
observation of the positive real part of the elastic forward scattering amplitude 
of all hadronic processes at high energies. This is a necessary condition for 
such process because in the commonly used formula for the refractivity index 
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its excess over 1 is proportional to this real part. However, later el I noticed 
that for such thin targets as nuclei the similar effect can appear due to small 
confinement length thus giving us a new tool for its estimate. If several gluons 
are emitted and each of them generates a mini-jet, the ring-like substructure 
will be observed in the target diagram. If the density of mini-jets within the 
xing is so high that they overlap, then they form a ridge pattern (or a wall 
0). If the -number of emitted gluons is not large, we'll see several jets (tower 
structure 13) correlated in their polar, but not in the azimuthal angle. Central 
collisions of nuclei are preferred for observation of such effects because of a 
large number of participating partons. If the number of correlated particles 
within the ring is large enough, it would result in spikes in the pseudorapid- 
ity distributions. However, the usual histogram method is not always good 
to verify these spikes because it may split a single spike into two bins thus 
diminishing its role. However, some hints to such structure can be found from 
these histograms. The wavelet analysis is well suited for this purpose because 
it clearly resolves the local properties of a pattern on the £vent-by-event basis 
by choosing the so-called Heisenberg windows (see, e.g., El). 

Another example is provided by the elliptic flow i.e. the azimuthal asym- 
metry Jn individual events. It may be related to a collective classical sling- 
effect lI of the rotation oil-colliding nuclei after peripheral collisions initiated 
by the directed pressure t3 at some impact parameter or to some jetty struc- 
tures. It can tell us about the equation of state of the hadronic matter by the 
shape of the created squeezed states. By passing, let me say that we also 
observed elliptic flow patterns ("cucumber") corresponding to the large value 
of the second Fourier coefficient, and the "three leaves flower" pattern with 
large third coefficient but I do not discuss them here. 

The event-by-event analysis of patterns in experimental and Monte Carlo 
events becomes especially important with the advent of RHIC and LHC. The 
4ir acceptance of detectors will be crucial for the event-by-event analysis . 

When individual events are imaged visually, the human eye has a tendency 
to observe different kinds of intricate patterns with dense clusters (spikes) and 
rarefied voids. However, the observed effects are often dominated by statistical 
fluctuations. The method of factorial moments was proposed til to remove the 
statistical background in a global analysis and it shows fractal properties even 
in the event-by-event approach- nn 

Both in cosmic rays EJO'EjtafiJ and at accelerators E3lia some ring-like 
structures were observed in individual events. Especially astonishing was the 
famous NA22 event of pion- nucleoli interaction with the pseudorapidity spike 
60 times exceeding the average level. 

Some event- by- event attempts (see, e.g., Eao) to treat nuclei collisions by 
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methods of the traditional inclusive and correlation measures revealed that 
there are some jetty structures in individual events which lead to spikes in 
the angular (pseudorapidity) distribution. Various Monte Carlo simulations 
of the process were compared to the data and-faUed to describe this jettyness 
in its full strength. More detailed analysis EJEj of large statistics data on 
hadron-hadron interactions (unfortunately, however, for rather low multiplic- 
ity) with dense groups of particles well separated from other particles in an 
event showed some "anomaly" in the angular distribution of these groups. 
It was found that the centers of these groups prefer to be positioned at a 
definite polar angle. This feature favors the above interpretation in terms of 
Cherenkov gluons. 

The wavelet analysis reveals the local properties of any pattern in an 
individual event at various scales and, moreover, avoids smooth polynomial 
trends and underlines the fluctuation patterns. By choosing the strongest fluc- 
tuations, one hopes to get those dynamical ones which exceed the statistical 
component. 

First attempts tn use wavelet analysis in multiparticle production go back 
to P. Carruthers cJO'Eil who used wavelets for diagonalisation of covariance 
matrices of some simplified cascade models. The proposals of correlation 
studies with the help of wavelets were promoted uE-L and used, in particular, 
for looking for the disoriented chiral condensate rW% The wavelet transform 
of the pseudorapidity spectra of JACEE events was done in c3. 

The wavelet analysis of patterns in central Pb-Pb collisions at 158AGeV 
was first done in 123'til. Five events with the highest registered multiplici- 
ties from 1034 to 1221 charged particles were chosen from 150 events of the 
emulsion chamber experiment EMU15 at CERN by the group from Lebedev 
Physical Institute. 

The target diagrams of secondary particles distributions for these events 
are shown in Fig.l, where the radial distance from the center measures the 
polar angle 9, and the azimuthal angle <fi is counted around the center. One 
can sum over the azimuthal angle and plot the corresponding pseudorapid- 
ity i 7 l = — logtan^/2) distributions shown in Fig. 2. The pronounced peaks 
(77-spikes) strongly exceeding expected statistical fluctuations are seen in in- 
dividual events. This inhomogeneity in pseudorapidity can arise either due 
to a very strong jet i.e. a large group (tower) of particles close both in polar 
and azimuthal angles or due to a ring-like (ridge) structure when several jets 
with smaller number of particles in each of them have similar polar angle but 
differ ipjtheir azimuthal angles. 

In wavelets were first used to analyze two-dimensional patterns of fluc- 
tuations of the event 19 in Fig.l. The results of the one-dimensional analysis 
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E3 of the two-dimensional target diagram are presented. To proceed in this 
way, the whole azimuthal region was divided into 24 sectors. The pseudora- 
pidity distributions in each of them were separately analyzed. Neighboring 
sectors were connected afterwards. Both jet and ring-like structures are found 
from the values of squared wavelet coefficients as seen from Fig. 3 E3. At small 
scales a, the wavelet analysis reveals individual particles. At larger scales, 
the clusters or jets of particles are resolved. Finally at ever larger scale one 
notices the ring-like structure around the center of the target diagram which 
penetrates from one azimuthal sector to another at nearby values of the polar 
angle (pseudorapidity) , thus forming an elliptic ridge. This structure approx- 
imately corresponds to the peak in the pseudorapidity distribution (for more 
detail, see Ej). Let us note, that it is not easy to notice in the target diagram 
of this event any increase of density within the ring just by eye because of the 
specific properties of the 9 — <fi plot where the density of particles decreases 
fast toward the external region of large polar angles. 

To reveal these patterns in more detail one should perform the two- 
dimensional local analysis. It is strongly desirable to get rid of such drawback 
of the histogram method as fixed positions of bins that gives rise to splitting 
of a jet into pieces contained in two or more bins. This chance is provided 
by the wavelet transform of particle densities on the two-dimensional plot. 
Wavelets choose automatically the size and shapes of bins depending on par- 
ticle densities at a given position. 

In principle, the wavelet coefficients Wj ly k 1 ,j 2 ,k 2 of the two-dimensional 
function f{6, 4>) are found from the formula 



Here Oi, <pi are the polar and azimuthal angles of particles produced, f(6, (f>) — 
5(8 — 6i)5(4> — (pi) with a sum over all particles i in a given event, (fci, £2) 
denote the locations and (ji,j2) the scales analyzed. The function ip is the 
analyzing wavelet. The higher the density fluctuations of particles in a given 
region, the larger are the corresponding wavelet coefficients. 

In practice, the discrete wavelets obtained from the tensor product of 
two multiresolution analyses of standard one-dimensional Daubechies 8-tap 
wavelets were used. Then the corresponding^s, sd and dd coefficients in the 
two-dimensional matrix were calculated (see a) . The common scale j% — fa = 
j was used. 

As stressed above, the ring-like structure should be a collective effect 
involving many particles and large scales. Therefore, to get rid of the low- 
scale background due to individual particles and analyze their clusterization 
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(1) 



properties, we have chosen the scales j > 5 where both single jets and those 
clustered in ring- like structures can be revealed as is seen from Fig. 3. There- 
fore all coefficients with j < 6 are put equal to 0. The wavelet coefficients for 
any j from the interval 6 < j < 10 are now presented as functions of polar and 
azimuthal angles in a form of the two-dimensional landscape-like surface over 
this plane i.e. over the target diagram. Their inverse wavelet transform allows 
to get modified target diagrams of analyzed events with large-scale structure 
left only. Higher fluctuations of particle density inside large-scale formations 
and, consequently, larger wavelet coefficients correspond to darker regions on 
this modified target diagram shown in Fig. 4. Here we demonstrate two events 
(numbered 3 and 6) from those five shown in Figs.l and 2. They display 
both jet and ring-like structures which are different in different events. To 
discard the methodical cut-off at r\ « 1.6 — 1.8 the region of 77 > 1.8 was only 
considered. 

Even though the statistics is very low, it was attempted to plot the pseu- 
dorapidity distribution of the maxima of wavelet coefficients with the hope to 
see if it reveals the peculiaritiesxibserved in high statistics but low multiplicity 
hadron-hadron experiments Ej'E3. In Fig. 5, the number of highest maxima 
of wavelet coefficients exceeding the threshold value Wj k > 2 • 1CP 3 is plot- 
ted as a function of their pseudorapidity for all five events considered. It is 
quite peculiar that positions of the maxima are discrete. They are positioned 
quite symmetrically about the value i] w 2.9 corresponding to 90° in cms as it 
should be for two Pb nuclei colliding. Difference of heights is within the error 
bars. More interesting, they do not fill in this centraLcegion but are rather 



For comparison, there were generated 100 central Pb-Pb interactions with 
energy 158AGeV according to Fritiof model and the same number of events 
according to the random model describing the inclusive rapidity distribution 
shape. The fluctuations in these events are much smaller than in experimental 
ones and do not show any ring-like structure. 

Thus I conclude that even on the qualitative level there is the no- 
ticeable difference between experimental and simulated events with larger 
and somewhat ordered fluctuations in the former ones. My aim here 
is to show the applicability and power of the two-dimensional wavelet 
analysis , the qualitative features and differences leaving aside quantita- 
tive characteristics till higher statistics of high multiplicity AA-events be- 
comes available. In particular, the special automatic complex for emul- 
sion processing with high space resolution in Lebedev Physical Institute 
(www.lebedev.ru/structure/pavicom/indcx.htm) is coming into operation, 
and it will enlarge the statistics of central Pb-Pb collisions quite soon. STAR 
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Collaboration of RHIC should also give soon some data at ever higher energy. 
Wavelets provide a powerful tool for event-by-event analysis of fluctuation 
patterns i n such collisions. i 
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The target diagrams of five events of central Pb-Pb collisions at energy 
158 GeV/nucleon obtained by EMU-15 collaboration. 
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The pseudorapidity distributions of particles in five events shown in Fig.l. 
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Fig. 3 The wavelet coefficients for the event 19 analyzed in El. Dark regions 
correspond to large values of the coefficients. Four of 24 sectors are shown. 
Rapidities are along x-axis, the scales increase down the vertical axis. 



radmOO: submitted to World Scientific on February 1, 2008 8 



Fig. 4 The modified large-scale target diagrams of two events (3 and 6). 
Darker regions correspond to larger particle density fluctuations. 
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Fig. 5 The pseudorapidity distribution of the maxima of wavelet coefficients. 
The irregularity in the maxima positions, the empty voids between 
them and absence of peaks at r] w 2.9 are noticed. 
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